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Abstract. The communications antenna on a spin stabilized spacecraft is seldom located on the spin axis, 
hence, the antenna is in motion relative to the center of mass of the spacecraft. The Doppler shift observed 
at the ground or space relay communications receivers will include oscillations whose frequency and 
amplitude are functions of the motion of the antenna and the attitude of the spacecraft relative to the line of 
sight (LOS). This functional dependence creates the possibility of estimating attitude parameters from 
Doppler measurements. This paper presents mathematical models of Doppler oscillations from spinning 
spacecraft, including the effects of nutation. Algorithms for estimating spin rate, attitude and nutation angle 
are described. Results of analysis of Doppler tracking of GOES-8 and WIND are also discussed. 

Introduction. Doppler tracking measurements are usually treated as equivalent to measurements of the 
range-rate of the spacecraft center of mass. However, Doppler measurements actually measure the range 
rate of the antenna, which during attitude maneuvers or spin stabilization has its own motion relative to the 
center of mass. A simple model for a spinning spacecraft is to assume that the angular momentum and 
angular velocity vectors are collinear. Then, unless the antenna location is on the spin axis, the antenna will 
have uniform circular motion relative to the spacecraft center of mass, so that the range rate will vary 
sinusoidally. The frequency of this motion is the spin rate; the amplitude of the motion is a function of the 
angle between the LOS direction and the spin axis. Hence, Doppler data can be used to estimate the 
orientation of the spin axis. There are two major phenomena that may invalidate this simple model. One is 
nutation, where the angular velocity vector is rotating about the (inertially fixed) angular momentum vector. 
The other is multipath, where the propagation path from the spacecraft antenna to the tracker includes 
reflections from other parts of the spacecraft, parts which have different Doppler shifts than does the 
antenna. In addition, a spacecraft may have more than a single antenna. Multipath will not be discussed in 
this paper. 

Apparent amplitude of Doppler oscillations. The amplitude of the oscillations in Doppler tracking 
caused by satellite rotation is a function of the antenna location, the rotation rate, and the aspect angles 
relative to the uplink and dowmlink trackers. For the purpose of deriving this functional dependence, 
effects of nutation will be ignored, hence, the spin axis is assumed to be coincident with the angular 
momentum vector, which is fixed in inertial coordinates. The inertial coordinate system will have the spin 
axis as the z-axis, and the x- and y-axes orthogonal. The following parameters will be used, 

s - -unit vector along spin axis 

u - -unit vector along LOS from uplink tracker to spacecraft in inertial coordinates 
u = iw, + j u 2 + kw 3 

d - -unit vector along LOS from downlink tracker to spacecraft in inertial coordinates 
d = \d x + j d 2 + k d y 

r A - -distance from spin axis to antenna 
s--unit vector along spin axis 

r(f) - -position vector of antenna relative to spacecraft center of mass in inertial coordinates 
r(0 = ir, (r) + jr 2 (r) + kr 3 


*This work was supported by the National Aeronautics and Space Administration (NASA)/ Goddard Space 
Flight Center, Greenbelt, Maryland under contract NAS 5-31000. 
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v(f) - -velocity vector of antenna relative to spacecraft center of mass in inertial coordinates 
radians 

( o s - -spin rate in 

sec 

The position and velocity of the antenna as functions of time can be written, 
r(f) = i r A cos(co 5 /) + \r A sin((0 A f) + k/* 3 

d i v 

v(f) = —?(/) = r A (& s \—ism((O s t ) + jcosfco/)) 
at 

The projections of the LOS vectors on the antenna velocity are, 

u • v(/) = (u - (u • s)s) •\(t) = r A co 5 {-u ] sin (D a f + u 2 cosco/) 

d* v(0 = (d-(d*s)s)* v(/) = r A (x) s (-d ] sin a )J + d 2 costa/) 

Where c is the speed of light in the same units as range rate and fj is the transmit frequency in Hz , define 



K is the turnaround ratio of the spacecraft transponder: K = 240/221 for NASA S-band trackers. The range 
rate oscillation and associated frequency deviation observed at the downlink tracker are, 

R A (t) = r^co j (~(w, + d,)sinco ; + (w ; + rf 2 )cos<o/) 

/.(/) = R,(t)A f 

The expression for frequency deviation can be converted to amplitude-phase form, 

/ A (0 = ^cos(co,r + <t>) 

The amplitude of the frequency oscillation is, 


r A^A f v( 

“waI 

u, +d ,) 2 +(«, +d ,) 2 
u + d-((u + d)*s)s| 

= 

/flflf+HI +2 (u*d)-((* i + d)*s)' 

= r A®Af\ 

/2 + 2(u«d)-((u + d)«s) 


For 2-way tracking, where the transmit and receive trackers are the same, the amplitude reduces to 
A = 2r A (0 J A / 7l-(u*s) 2 

The 1-way amplitude would be half the 2-way. Note that there is a 180 deg ambiguity; if the spin vector is 
replaced by it's negative, i.e., the direction of spin is reversed or the spacecraft is flipped over, the apparent 
Doppler oscillation amplitude will be the same. 

Estimating attitude. Given that the spin rate is known, all the parameters in the above equations are known 
except the spin axis unit vector. The absolute value of the projection of the spin axis unit vector onto a 
known vector can be solved for as follows: 
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(2 - way) |u • s| = 


(1 - way) 


(3 - way) 

The last equation can be rewritten in terms of unit vectors by normalizing, 

(u + d) „ 

(3 - way) t; Hr • s 

u + d 

From this point on, it will be assumed that the Doppler tracking observations have been reduced to a set of 
projections and unit vectors, 

o*=|*t #s |; 0<a k <\ 

Each projection is equivalent to an angle, 

S* = COS' 1 (a,); 0<^<1 

When many observations (projections and unit vectors) are available, a least squares solution can be 
computed. Depending on geometry, several solutions are possible, so initiating the least squares algorithm 
properly is important. If the attitude is the only unknown, the right ascension (a) and declination (5) of the 
spin axis can be solved for directly by a non-linear least squares routine, 

cosScosoc 
cosSsina 
sin5 

If another parameter such as antenna radius from spin axis is uncerL.’ may also be solved for. The 
MATLAB™ Optimization Toolbox installed on FDF LANs contains built functions for minimize* 1 m 

Analysis of GOES-8 Doppler Tracking. GOES-8 was launched 13 April 1994. Doppler tracking o 
was received at FDF from DS46 (Canberra) 3-way with the Indian Ocean Station (10S or SEYS) teleme y 
system as the transmit site. For data taken from approximately 0700Z to 0740Z, the spacecraft w f as spin 
stabilized. Figure 1 show’s the Doppler oscillations from valid data obtained during this time. Shown by a 
dotted line is the estimate of amplitude derived from these oscillations. The measured rotation rate is 
0.0182 Hz, while the distance of the antenna from the spin axis is assumed to be 3.3 meters. Figure 2 
depicts contours of the surface generated by varying right ascension from 25 to 75 deg and declination from 
-50 to 0 deg and plotting the reciprocal of the resultant RMS error between the measured projection and the 
vector dot product. Further analysis show-ed that the fit could be improved if r^ was varied as well. The 

MATLAB routine LEASTSQ in the Optimization Toolbox was used to solve for right ascension, 
declination and antenna radius from the spin axis. The routine was run twice, initialized with each of the 
angle pairs found above and with r^ initialized at 3.3 meters. The solutions were, 




r »/( 2 + 2 ( fi # d ) 
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RMS error 

0.042 

0.044 

Right ascension (deg) 

62.2 

59.7 

Declination (deg) 

- 24.7 

- 16.7 

Antenna radius from spin axis (m) 

2.79 

2.76 


GCES-t DoDpte*- osollaOoos and envelope 



Figure 1. GOES-8 Doppler Tracking. 


Since the observations were not exactly co-planar, the solution with the smaller RMS error would be the 
more likely candidate for the true attitude. The planned attitude after separation was a = 61.7 deg and 5 
= -25.0 deg, according to reference 1. Figure 3 shows the fit obtained by comparing the projection of the 
solved for spin attitude on the unit vectors with the projections measured from tracking data. 


contours of reciprocal error surface 
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corTputed pfOf®cl>on arxl measured pcojecttor 



Figure 3. Measured versus computed projections. 


Spinning spacecraft attitude dynamics (nutation). Analysis of the dynamics of a rigid spinning 
spacecraft in order to predict the effects of such motion on Doppler tracking data uses the following 
parameters: 

1. The spacecraft major moments of inertia, in Kg-m^. The nominal axis of rotation is assumed to 
be the z-axis, so the values are represented as a vector [lx Iy Iz], where the labels x, y and z are 
assigned so that the entries are in monotonic (lowest to highest or highest to lowest) order. The 
choice of the axis with the intermediate moment of inertia as the spin axis is unstable and will not 
be considered. 

2. The location of the spacecraft antenna in meters, as a vector [x^ y^ z^] in body coordinates. 
The assignment of x^, y^ and z^ should be consistent with the moments of inertia. 


3. The angle (3 of the LOS vector with the angular momentum vector. 

4. The angular momentum magnitude L. 

5. The initial nutation angle 0 q, which is assumed to be small. 


Reference 2 is the source of the basic attitude dynamics relationships in the following. The convenient 
inertial reference frame has the angular momentum vector along the z-axis, with the x and y orthogonal. 
The time reference is chosen so that at t=0, the angular velocity vector is in the body x-z plane, so the 
angular velocity vector in body coordinates is, 


®*0 =i/,w.o + j/. 0+ k*t>) 3 o =i 


Lsin0 o u Lc os0 o 

■ I K L 

lx Iz 


CO 


10 


lx ' 


“30 = 


L_ 

Iz 
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The kinetic energy E is. 


lx cojp + Iztolo 1 f L 2 sin 2 0 O L 2 cos 2 0 C 


The "parameter'’ m is given by. 


(!x-Iy )(L 2 -2 IzE) V lx h 

(■ lz-ly)(L 2 -2IxE ) /f fw2 f, ,fsin 2 e„ cos 2 9 0 

( Iz- h)L \-lx - + - 

' l lx h 


v /a* 


Iz(Iy-Ix) 


Uz-fy) I 


(fe-/jr)cos A 9 0 ] /x(/:~/v) 


Define the “normalized" parameter, 


Iz(Iy-Ix) , a2 
= ~ m / 0 n • 


/a(/z - /v) 


The body nutation rate is, 


CO = ± 
/> “1 


j (/:- /y)(L 2 -2IxE) , l(h- ly)L 2 ((h- Ix)cos 2 9 f 
' Ixlyh ~ V /.v/v/r 


==+J L \ Uz-ly)(lz-Ix) 


lxlylz 

\{lz- Iy){lz- lx) 


where the "+" is used if lx > Iy > Iz and if lx < Iy < Iz. 

Finally, the amplitude of the y body component of angular velocity is. 


C0 20 =, 


L 2 -2 IzE (Ix-h) sin 2 0. 


(ix-h) 


/>’(/>• - Iz) 


Ixlyily - ft) 


y Ixly(ly-Iz) 


_ L9 » IEUzML - L9 n I IxJ y - Aft - /xft + Zxft _ L8 ,, | 7-(/x-/v) 

A’ \lx{Iy-Iz) Iy \ /jr(/v-ft) ~ fy \ + Ix(ly-lz) 


= T 1 V e o + m = 7 1 0 oV 1 + m n 

/y /v 


The angular velocity vector as a function of t (sec) is given by, 


a), = co 10 cn(co p rlm) 
to 2 = -co 20 sn(co /7 /Im) 
co 3 = a) 30 dn(to rl™) 
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where sn, cn and dn are the Jacobian elliptic functions. When m « 1, approximations for the Jacobian 
elliptic functions are, 


sn(x\m) = sin(x) 
cn{x\m) = cos(a') 

dn{x\m) = 1- — sin 2 (;t) = l - — + — cos(2a*) 

2 4 4 

So the angular velocity vector in body coordinates can be written, 


to, = co l0 cn(w r fl»i)* a> l0 cos(ay) 
to 2 = -<o 2i) sn(<& p t\m) ~ -o) 2u sin(co p /) 


co 3 = (o 30 dn(co rim) ~ co 3( 


m m 
4 4 



The instantaneous velocity of the antenna in body coordinates is found by crossing the angular velocity 
vector with the antenna position vector. The antenna velocity is then rotated to inertial coordinates and 
projected on the LOS vector to give the range-rate of the antenna as a function of spinning. 

The Euler rotation angles for angles 0 and \\f can be found by solving. 




"o' 


sin Bsin \y 

^b2 

II 

s 

CD 

k 

0 

= L 

sin 0cos \\f 



L_ 


cos0 




Lsin 0sin \\f 



Ix 

03, 


Lsin0cos\]/ 





Iy 

©3. 


L cosG 



Iz 


So the minimum and maximum values of nutation angle are, 


0 

0 


mm 


max 


= 0 (. 

— ^0 iti ~ 0 


This implies, 

6 L ^ e o + m = 1 ! ffl _ b'Uz - Ix ) 

0o e 2 0 2 lx{h-Iy) 

which agrees with equation 16-114 in reference 2. 


The angular velocity vector in body coordinates can now be written as, 
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f>, =— e 0 cn((O rim) 
lx 

U>2 m 

ly 

L 

o^- 


Sine and cosine of the Eulerian angle i| / can be found as functions of time by solving. 


tan\j/ = 


lx CO, 


sinxjr = 


/ A'CO, 


i -;cos\|/ = - 

^ VC0 2 ^Jlx cOj + Iy 2 co 2 2 - yjlx 2 o>* + Iy 2 a > 2 2 


ly(0 2 


/_r 2 oo ^ + /v 2 co 5 = L9 0 ^c^(co p rlm) 2 + (l +m 0 )^n(co p Hm) 2 
= L0 O ^/l + m 0 5/7(co p /lm) 2 = L0 


sin \j/ = 


cosy - 


/jCCOj 


^jlx 2 (o 2 + /y : to 2 


/jt( — )c«(co rim) 

lx ) 0 O . , . 

= — c/Kco/lm) 

0 


/yco 2 




L0 


-*£"• J sn<ra / 1 "’) _ e 


-Jb? cof + /y : co 


z.0 


s/i(co/lm) 

0 


A differentia] equation relating $ to 0 and y is, 


4> = (cOj sin y + co 2 cosy)/ sin 0 = L 


/ • 2 2 

' sin \|/ cos \j/ 

/* /7 


« L 


2 ^ 


( 


J_ 0Q (lx- Iz)(!x~ Iy)sn(<£> p t\m) 
lx Q 2 Ix 2 (Iy-Iz) 

1 0 2 {lx - Iz){Ix- Iy)sn(co p t\m) 2 

lx ( 0 ^ + msn((O p t\m) 2 )lx 2 (ly - Iz) 


L L(lx- 7v)cos 2 y 
lx Ixly 


To get the frequency of oscillation of terms involving cost)) and sin<|>, the average is needed. 


1 

Ix + \ 

' Q\{Ix - lz)(lx - ly ) > 


■ sn((ti p t\m) 2 dt 

l Ix 2 (Iy-Iz) j 

2K } 

K (e 2 0 +msn(d) p t\m) 2 ) ^ 


where K is the half period of the Jacobian elliptic functions (ref. 3, section 16). The integral may be 
approximated by substituting sin for sn since m is small. 
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1 c sn 

2K_Uel + 


sn((0 t\m) 2 dt 


<Q P 7 sin 2 (o i p t)dt 

TT J ( A 2 


= L J 


+msn((& p t\*n) 1 ) n L (Op + m - mcos 2 (co„r)) 
sin 2 (x)dx 




+ m - mcos 


1 

, run 

>)”« 

l i +m ) 


(integration formula 262 in ref. 4) 
Substituting for the integral. 


< 4> >= L 


_ 1 _ 

Zr 




9S(/jr-fe)(/r-/.v) ) 1 


/jC 2 (/> - /Z) 


m 


1 - 


6 o 


' 6 q + m 


Jx 


m 


= L 


= L! 


Av 


( /x - /')( /.v - !y) 1 8 ' 

l Ix 2 (Iy-Iz) j 

Ux-!z)Ux-ly) ' 
lx 2 (ly- Iz) 

( 


1 - 


V 


1 + ml 9 ( : , j 


Ixilz-Iy) 


( 


lx 


J_ + [ (lx- Iz) 


lxlz 


1 - 


V 


lz(Iy-l. x) 
W 


1 - 


V 


Ix(lz-ly) 

ly(Iz-lx) 


J , 


Ix(Iz-ly) 

ly(Iz-lx) 




+ L Mt-lxWz-ly) 
Iz V Wv 


L (Iz- Ix)(Iz- ly) 


Ixly 


; /a* > /z 


> = -®r 


Thus, the signals cost)) and sintj) oscillate at an average frequency approximately equal to the difference of 
the spin rate and the nutation frequency. 


Estimation. With the preceding derivation, it is now possible to write the component of measured 
range-rate caused by antenna motion as a function of known parameters, signals with known frequencies 
and amplitudes and parameters to be estimated. The parameters are: 


parameter 

Name 

known or to be estimated 

L 

angular momentum magnitude 

either 

0 

nutation angle 

to be estimated 

P 

LOS angle 

either 

lx, Iy, Iz 

moments of inertia 

known 

mo 

normalized parameter of elliptic 
function 

known 

X A* >'A> Z A 

body coordinates of antenna 

known 
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v (0 = -j 

dt 



1 

'c * 

I 


1 

1 


d_ 

dt 

d_ 

dt 


(cosy cos<{) -siny sin 4 ) )a- 4 +(~siny cos(|)-cosy sin<{)) v A + 0z A cos § 

x 

0 sinyx A + 0 cosy v A + z A 

cos(<|) +y )x A + -sin(<j)+y )y A + 0cos<j)z A 
x 

0 sinyx A +0 cosyy A + z A 


The second component will not be used and so is not calculated, since the LOS vector in inertial 
coordinates is assumed to be. 


LOS = 


sin (3 
0 

[_COSp 


The measured range-rate will be the dot product of the LOS vector and the inertial velocity, hence the 
incremental range-rate due to antenna motion that will be measured by the tracker is the sum of two 
components. The two components will be evaluated separately. 

If the LOS direction w-ere parallel to the angular momentum vector, the range-rate due to antenna motion is, 
u'here <I> is an unknown phase angle, 

R v *-j-(0sin\| fx A + 0 cos vj/y^ + z A ) 

at 

* 9 ° Cn (°V lm ) “ >’>( 0 max ^ + Z A ) 

= (-*,e 0 ™(cvlm) - y A Q ^cnica p t\m))a> p dn((£> p t\m) 

This will be called the "vertical component" of the antenna motion. It can be rewritten as, 

R v = (-.v„9 0 .m(a) / ,fl m) - y A % /l + m (t B 0 cn(oy p t\ m)^ r dn((a r tl m) 

= 00 + U + '«o ^ COS(C 0 / + O) 


So the vertical component is a signal oscillating at the nutation rate and with amplitude proportional to the 
minimum nutation angle. 

The horizontal component is a bit more complicated. 
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R h ~ — [(cosy cos $ - sin V|/ sin 0 ) + (- sin y cos 4 > - cosy sin <>) y A + z^Ocos^] 
dt 


-0 0 t 

— — 5fi(0) p /lm)cos<t) — - cn ( to p i\m) sin $ | + 
0 0 




0 0 

— 1 — — -c/i(co p t\m) cos $ | + z A 0cos<j 

0 0 


0 - 


(x A {-j\ + m 0 sn(w p t\m)cosQ>-cn((O p t\m)s]n$)+ y A {J\ +m 0 sn(Vi p t\m)s\nQ-cn(<a p t\m)cos$)) + z A Qcos$ 


The ratio 0 O / 0 has essentially no dependence on 0, so the nutation component of antenna motion is, where 
<J>' and d>" are unknown phase angles, 


= 4 ^ 9COS(}) ) 

dt 

= z A ( 0 cos <() — 0<b s i n <t> ) = ’. 4 (ecos<j)- 0(co, o - co ;i )sin6) 

= -zA + -° } ,.)s«n((0) ;(0 -co ,,)/ + <*>') + 

<-A d o ,n o .1(0) +0 ) )sin(((o v , + to )/ + O') -(co, 0 -3co / ,)sin((o>- 0 - 3co,, )r + 0')] 
2/4 + 2 m 0 

So ihe horizontal nutation component is a signal oscillating at a frequency equal to the difference of the spin 
rate and the nutation rate, together with sidebands at intervals of ±2otp, all with amplitudes proportional to 

the minimum nutation angle. 

WIND Nutation Estimation. During initial phases of FDF support of the WIND mission, Doppler 
tracking data was used to monitor the spin rate and nutation angle of the spacecraft. The purpose was to 
provide a back-up to on-board attitude sensor telemetry, with the emphasis on measurement of nutation 
angle immediately prior to the Lanyard Boom deployment. The WIND spacecraft was launched on 1 
November 1994. WIND is spin stabilized with a nominal spin rate of 20RPM. WIND is equipped with a 
nutation damper with a nominal time constant of 1400 seconds. The spacecraft is cylindrical, with the z- 
(spin) axis being the axis of the cylinder. The Medium Gain Antenna (MGA) is extended from the negative 
z face of the cylinder. Figure 4 show’s the WIND spacecraft configuration. The WIND attitude is 
maintained w ith the spin axis parallel to the south ecliptic pole, effectively perpendicular to the LOS. 
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The nutation estimation results for the first (of three) thruster flush burns will be presented. During the 
maneuver, the WIND moments of inertia were 741, 678 and 912 Kg-m 2 . The MGA body coordinates are 
(0.14,-1.15,-2.68) meters. For WIND’S attitude, according to the analysis above, the Doppler signal will 
contain a component at a frequency approximately equal to 1.3 times the spin rate, with an amplitude 
proportional to the minimum nutation angle (the nutation angle oscillates between minimum and maximum 
values since WIND is not axially symmetric). For the assumed WIND configuration, the minimum nutation 
angle in degrees can be estimated by finding the amplitude of this component in meters/sec, and dividing it 
by (0.007SR), where SR is the spin rate in RPM. 



detrended Doppler observations 









The topmost plot in figure 5 shows WIND range-rate for the period of the thruster flush maneuver. The 
middle plot shows the detrended data, used as the input to a Discrete Fourier Transform (DFT) algorithm. 
The lowest plot is the recovered nutation signal obtained by performing the inverse DFT on only those 
frequency components near the expected location of the nutation signal. Since some frequency components 
of the desired signal are lost, and some of other signals are included, there are large oscillations at the 
beginning and end of the recovered signal; these are known as the Gibb’s Phenomenon (ref. 5, pp. 73-75). 
The nutation signal was “envelope detected” and scaled by 1/(.007SR) to obtain the minimum nutation 
angle (0 o ) as a function of time. The comparison is with Sun Sensor Assembly #2 on WIND (listed as 

SSA2 in figure 4.). SSA2 is 30 degrees away from the X-axis, so the largest nutation angle it sees is 
between 0 () and 0 MAX . The Doppler nutation signal was further scaled by , 


|i+Sl 

m Q 

— cos 

2 ( 

' 30 k y 

1 2 

2 

L \ 

,180 


;/?7 0 = 0.49 


to put the SSA2 and Doppler signals on the same reference. Figure 6 show's the comparison of the 
detrended SSA2 angles and the scaled Doppler nutation signal amplitude. 


SSA2 and Doppler derived nutation ang!es--WIND 04 NOV 1995 



Figure 6. Comparison of Sun Sensor 2 and Doppler Nutation. 


Conclusions. For spinning spacecraft, processing of Doppler tracking data for attitude parameters can be a 
useful adjunct to on-board attitude sensors. For “smaller, cheaper” missions it could possibly replace 
on-board sensors. In order to be useful for real time attitude estimation, further algorithm development 
should be done; for example, a digital phase-locked loop (DPLL) implementation for tracking spin rate and 
demodulating the nutation signal. 
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